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ABSTRACT 



2. BUILDING THE RESPONSE MATRICES 
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We have used the IBIS Simulator to produce prelim- 
inary response matrices for the ISGRI and PICsIT 
detectors in order to help understanding their scien- 
tific performances before the calibration results are 
available. The derived matrices, in a format compat- 
ible with the XSPEC spectral analysis package, have 
been tested by fitting simple models and then used 
to analyze simulations of astrophysical sources with 
more complex spectra. 

Key words: IBIS; Monte Carlo Simulations; Instru- 
mentation. 



1. THE SIMULATOR 



The IBIS Simulator (developed in collaboration 
by IFC/CNR Milan, Southampton University and 
ISDC) can perform the following steps. It first gen- 
erates a model of the gamma-ray sky with the pos- 
sibility of defining position, intensity and spectral 
properties of the celestial sources. Then different ob- 
servation strategies (pointing directions, observation 
length, dithering patterns, etc.) can be defined. Af- 
ter these preliminary steps the interaction of all the 
source photons with the active and passive materials 
of the instruments are simulated. An estimate of the 
background is also added, based on the results of the 
INTEGRAL Mass Model (Lei et al. 1999), which is 
a detailed Monte Carlo simulation that takes into ac- 
count the interaction of the particle and photon flux 
with the whole INTEGRAL spacecraft. The Sim- 
ulator finally produces a data set in a format fully 
compatible with the ISDC analysis system. For the 
analysis of the output data we have developed some 
prototype programs, which perform image deconvo- 
lution and spectral extraction. 



To build the response matrices of ISGRI and PICsIT 
(Ubertini et al. 1999), two sets of simulations of on- 
axis monochromatic sources have been performed at 
different energies (20 values for ISGRI between 15 
keV and 600 keV and 18 for PICsIT between 100 
keV and 5 MeV). An example of the resulting data 
for some sample energies can be seen in Figures [l] 
and@. 



ISGRI calibration simulations 
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Figure 1. Simulated ISGRI response to on-axis 
mono- energetic sources of different energies. The 
spectra are binned with different factors. 

The simulations clearly show the photopeak, the 
Compton edge and, in the case of PICsIT the 
backscattering peak (see Figure ||). These three 
features have been fitted separately with analytical 
functions. The photopeak required a Gaussian func- 
tion, while the other two required a Gaussian plus a 
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PICsIT calibration simulations 
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Figure 4- Interpolation of the photopeak parameters 
of the ISGRI response. 



Figure 2. Simulated PICsIT response to on-axis 
mono- energetic sources of different energies. 

quadratic function. 

The resulting set of parameters as a function of the 
energy have been then interpolated over the entire 
range of energies of the two instruments (see Figure 
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PICsIT response to 1.5 fJeV photons 
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Figure 3. Simulated PICsIT response to a mono- 
chromatic source of 1.5 MeV photons. 

Finally these results have been used to create two 



FITS files (for each detector): an RMF file contain- 
ing the normalized matrix and an ARF file which 
contains the information about the efficiency of 
the instruments. These two files have been created 
in such a way that they are compatible with the 
XSPEC data analysis package (OGIP standard ver. 
1992a). 

The matrices have been tested with different input 
spectra. An example is the following: a ^lOOmCrab, 
on-axis source with a power law spectrum observed 
for 1050 s (input parameters: photon index r=2, 
Axookev = 9xl0~ 5 phem" 2 s" 1 keV^ 1 ). The best fit 
with XSPEC yielded T = 1.987 ± 0.0483 , A 100 fc e y = 
(8.744 ± 0.6357) x 10~ 5 ph cm" 2 s" 1 keV" 1 , Reduced 
X 2 = 0.942394. The confidence contours of the fit are 
presented in Figure |5[ 



3. SCIENTIFIC APPLICATIONS 

3.1. Nova Muscae (GRS 1124-684) 

The X-ray transient GRS 1124-684 (Nova Muscae) 
was discovered in 1991 with the GRAN AT and 
GINGA satellites. Measurement of its mass function 
during quiescence established the presence of a black 
hole (Remillard et al. 1992). GRS 1124-684 was ob- 
served several times by the SIGMA telescope during 
its outburst (Goldwurm et al. 1992). An emission 
feature at ~500 keV was discovered during the last 
13 hours of the January 20th observation. 

We have simulated IBIS observations of this spectral 
feature, with exposure times ranging from 125 s to 
15,000 s, using as input parameters the values mea- 
sured by SIGMA. An example of a 2000 s simulation 
(see Table 0) is shown in Figure 0. The simulations 
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Figure 5. Confidence contours of the power law pa- 
rameters for a ~1000 s observation of a ~ WOmCrab 



Figure 7. IBIS line sensitivit y ob tained from the sim- 
ulations described in Section S.i. 



Nova Muscae - Exposure time: 2000 s 
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1E1740-3942 - Exposure time: 4000 s 
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Figure 6. 2000 s IBIS simulation of Nova Muscae. 

show that the ~500 keV feature can be appreciated 
(~3 a) for observations as short as 500 s. 

3.2. IBIS line sensitivity 

We have performed simulations with different expo- 
sure times and spectral parameters to determine the 



Figure 8. 4 ks IBIS simulation of 1E1740-2942. 



sensitivity to detect an emission line at ~ 500 keV 
over a continuum. Figure |7] shows the minimum ex- 
posure time that is required for a ^3 a detection 
of such a feature as a function of the line flux and 
width (a) . The continuum used is that of Nova Mus- 
cae (see Table ??). A feature like that observed in 
IE 1740-2942 during October 13-14 1990 with the 



Table 1. IBIS simulation of Nova Muscae - Exposure time: 2000 s. 





Input parameter 


Best fit result 


Photon Index 
Norm, at 1 keV 

Line energy 
Line width (a) 
Line flux 


2.42 

11.27 ph cm" 2 s" 1 keV" 1 
481 keV 
23 keV 
6xl(T 3 ph cm" 2 s" 1 


2.423 ± 0.0184 
10.68 ± 0.741 ph cm" 2 s" 1 keV" 1 

483.3 ± 2 keV 
20.66 ± 0.1199 keV 
(3.017 ± 0.242) xlO" 3 ph cm" 2 s" 1 



SIGMA telescope (Bouchet et al. 1991), for exam- 
ple, can already be detected with an exposure time 
of ^250 s. We have also performed some simulations 
of this object. The parameters, that we have used, 
were a Gaussian line (center = 480 keV, a = ^100 
keV, flux = 1.3xl0~ 2 ph cm~ 2 s _1 ) over a comp- 
tonization continuum (kT = 27 keV, r = 3.2). In 
Figure ^ we present a 4 ks simulated IBIS spectrum 
of 1E1740-2942. 



4. CONCLUSIONS 



Of course the matrices presented here, being ob- 
tained trough simulations and not by real calibration 
data, are intrinsically limited in their accuracy ow- 
ing to the simplifications inherent to the simulator 
programs. Nevertheless they can be used to analyze 
simulated data and to obtain an estimate of the ob- 
servation time required to achieve specific scientific 
objectives. 

Future work will include the following effects that 
are not yet implemented in the current version: 

1) Charge loss in ISGRI, caused by the different 
mobility of the charge/hole pairs in semiconductors; 

2) Pixel disuniformity and/or gain variations; 

3) Multiple interactions in PICsIT; 

4) Angular dependence of the spectral response; 

5) Background spatial disuniformity; 
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